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Abstract 

We study the alpha vacua of de Sitter space by considering the decay rate of the inflaton field 
coupled to a scalar field placed in an alpha vacuum. We find an alpha dependent Bose enhancement 
relative to the Bunch-Davies vacuum and, surprisingly, no non-renormalizable divergences. We also 
consider a modified alpha dependent time ordering prescription for the Feynman propagator and 
show that it leads to an alpha independent result. This result suggests that it may be possible to 
calculate in any alpha vacuum if we employ the appropriate causality preserving prescription. 
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I. INTRODUCTION 



Field theory in de Sitter space is an area of active interest as the inflationary epoch can 
be well approximated by de Sitter space. Further, de Sitter space may even be applicable to 
the current universe as it is known that the expansion is accelarating . Current models 
of inflation postulate a scalar field rolling to the bottom of its potential as the source of 
the energy density that drives the expansion and quintessence models attempt to do the 
same for the present accelaration observed in the universe. Further, in inflation it is the 
quantum fluctuations in the scalar field that seeds the perturbations responsible for the 
temperature anisotropy and large scale structure and observations are in agreement with 



the predictions of inflation [3|. The nature of the de Sitter invariant vacuum is tied to trans- 



Plankian effects |4] in inflation. For these purposes it is important to investigate the formal 
construction of field theory in de Sitter space. 

Field theory in de Sitter space and the question ofpossible vacua of free scalar field theory 



in de Sitter space have been investigated by 




One of the surprising results that 



emerged was a vacuum choice ambiguity that is absent in Minkowski space. They found that 
in de Sitter space there exists an entire family of de Sitter invariant vacua that can be used 
to construct a consistent free field theory. This vacuum ambiguity is discussed for general 
spacetimes by Long and Shore '01 using a wavefunctional formalism. If one insists that the 
vacuum reduce to the Minkowski vacuum in the limit that the expansion vanishes then a 
single member of the family is selected and is commonly denoted the Euclidean or Bunch- 
Davies vacuum 3]- The vacuum choice ambiguity is normally addressed in inflationary 
calculations by assuming that the vacuum is the distinguished Bunch-Davies vacuum and this 
leads to successful predictions of the temperature anisotropics. The question of what physics 
is responsible for the vacuum selection remains to be addressed and a number of papers 
address the use of the general vacuum with a cut-off to encode trans-Plankian effects llll. 

..t.™.t.g tKeo.es tKe ,ene.a. vaeua, howeve. a.e a,„ed U^U - be 



inconsistent usually due to the presence of intractable divergences. 

Generally the dynamics of inflation are studied treating the inflaton, the scalar field 
driving inflation, as a classical homogenous field and the gaussian approximation is applied 
to the quantum fluctuations about this configuration. As mentioned these fluctuations seed 
the perturbations in the universe. For a review of the quantum treatment of inflationary 
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dynamics see Of particular interest to this paper is the question of particle decay in de 
Sitter space and [3l investigates this questions assuming that the Bunch-Davies vacuum is 
the relevant one. The same problem is addressed in a somewhat more general context using 
dynamical renormalization in 



Similar studies 



16l . IitI i of interacting theories in general a-vacua, the de Sitter in- 



variant family, are plagued by non-renormalizable divergences that obscure any attempt to 
probe their properties. Although this is a disturbing feature it does not seem reasonable to 
simply discard them as they are perfe ctly valid vacua of the free theory. There have been 
attempts to address these divergences and in particular we take up a suggestion put 

forward by Collins and Holman in (2^ when we study the problem of particle decay. We 
find their prescription especially attractive as it yields a vacuum independent result. 

In this paper we will address the question of the slowly rolling inflaton decaying into light 
scalars in a general a-vacuum. We will present a brief introduction to the nature of the 
vacuum ambiguity of de Sitter space in section (jllj) and present the notation we will employ. 
In section pi (^ll we will present the formalism necessary for calculating in a time dependent 
background. The interaction lagrangian and the calculation of the decay rate in the Bunch- 
Davies vacuum follows in section pil|) . We then present the results in general a- vacua 
for conformally coupled and minimally coupled in sections (jlll A|) and piIB|) respectively. 
Adopting a prescription for addressing time ordering problems in the a- vacua, we repeat the 
calculation for the conformally coupled field in section pil Finally we end with some 
concluding remarks regarding a consistent treatment of the a- vacua and their significance. 

II. BACKGROUND 

The quantization of a free scalar theory in a fixed de Sitter background and the attendant 
vacuum choice ambiguity is described in 0, 0|. Here we sketch the procedure _and describe 
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our conventions. Additional details of the procedure can also be found in 

We will find it convenient to work both in cosmic time(|H) and conformal timeQ in which 
the metric is respectively given by 

drf — dx^ , , 
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where H, as usual, is the Hubble constant. The domains of the coordinates are t G [—00^ 00] 
and T] E [—00, 0] and they are simply related hj 1] = — ^7^-- 

A free massive scalar theory satisfies the Klein-Gordon equation 

[V' + m|] $ = 0. (3) 

Since de Sitter space has constant curvature, any coupling to curvature normally denoted 
by ^-R can be absorbed into the mass and hence we take m$ to denote an effective mass into 
which the coupling has been absorbed. 

As usual the field is expanded in creation and annhilation operators 



where 



[a^,a^^] = {2Txf5\p-q). (5) 
Equation Q implies that the mode functions satisfy the differential equation. 



Ukiv) = 0. (6) 



This equation has Bessel functions as solutions. In particular, if we demand that the modes 
match the Minkowski massive free scalar modes at short distances then, 

U^iv) = ^v"'H^^\kv\ (7) 

where, 



/9 rn? 

and Hu denotes the Hankel funciton of the second-kind. The normalization of the modes 
has been fixed by enforcing the canonical commutation relation for the field and its conjugate 
momentum which leads to a Wronskian condition for the mode functions. We denote the 
corresponding annhilation operators a'^ and the vacuum annhilated by them, commonly 
referred to as the Bunch-Davies vacuum, as \E). So far, quantization has proceeded much as 
it does in Minkowski space but we have had to make the reasonable but arbitrary assumption 
that the modes match Mikowski modes at short distances to fix the vacuum. 
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A. Alpha vacua 



When we relax the Minkowski matching condition we see that any norm preserving lin- 
ear combination of the mode function and its conjugate will satisfy (jH)) and preserve the 
commutation relations. The standard literature addresses this by considering a Bogolubov 
transformation of the creation-annhilation operators, 

«^ = ^4«f-e"*aV), (9) 



where < 0. The vacuum annhilated by all a| is denoted Rewriting (@]) in terms of 
we find the corresponding mode functions, 

Unv) = N^iU^{v) + e-U^{vT). (10) 

It remains to be shown that this set of vacua parametrized by a - the Bunch-Davies 
vacuum is recovered at the a —oo limit - are infact de Sitter invariant. This can be 
demonstrated by showing that the two-point functions can be written in terms of the de 
Sitter invariant distance between the two points. We will construct the two point functions 



in the next section and address the question briefly while details can be found in 
references therein. 



Ill and 



B. Green's functions 



To compute the diagrams we will eventually encounter in evaluating the decay rate we 
need to construct the propagators for the free theory. Since the metric we are considering 
is spatially fiat we can write 

(11) 

This function depends only on Z{x,x'), the de Sitter invariant distance, and hence the de 
Sitter invariance of the a-vacua. In a momentum representation we have 

Gkiv,v') = unv)u:ivr, (12) 

as the momentum space Wightman function. 
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From the Wightman function we can construct the Feynman propagator following the 
same procedure as in Minkowski space, 



C'ix^x') = i{a\T[^x)^x')]\a), (13) 



(27r)3 

GUV,V') = &{r]-v')G{v,v') + &{^'-v)Q{r]\v)- (14) 

Defined in this manner the Feynman propagator, as would be expected, satisfies the equation 
of motion with a point source, 

[Vl + ml]G-{x,x') = ^-^-^^^. (15) 



Although natural, this definition for the Feynman propagator is not the only one leading 
to a well-defined object, and further, it is perhaps not the most meaningful defintion due to 



su 



Dtleties associated with the a-vacua and causality. These issues are addressed in |21l. 123. 
25l | and an a dependent time ordering in the definition of the propagator is suggested. Here 
we will present the central idea and the details necessary for calculation. 

In de Sitter space there is a natural antipode map given in the conformal time coordina- 
tization by xa '■ {v^ ^) ~^ i^V^ We can then expand 

g'^ix.x') = iV4e?^(x,x')+e"+"*G£;(x',x) 

+e''GE{xA, x') + e'^^GEix', xa)) • (16) 

This result indicates that ordinary time-ordering will not impose the causality expected. 
Instead we can define a time-ordering prescription that ensures a causal definition for the 
Feynman propagator and it leads to a double source propagator, 

[V^ + ml] G"(x, x') = Aj-^-^j^^ + B f)''^;'''\ (17) 

\/-9[x) V-9{xa) 

where Aa and must be fixed by some prescription while maintaining 

lim Aa = l and hm = (18) 



a— >— oo o^— oo 



SO that the standard time-ordering is recovered in this limit. Doubly sourced propagators 
have been discussed in jisl, but without the modified time ordering they lead to non- 
renormalizable divergences. We choose the particularly simple scheme A^ = 1 and B^ = 0. 



This leads to a singly sourced propagator and in the path integral formulation all internal 
lines are naturally converted to the Bunch-Davies propagator. This conveniently takes care 
of the divergences that plague calculations in the general a-vacua. We should note however 
that the two-point functions are still a dependent. 



C. Schwinger-Keldysh formalism 

Since de Sitter space presents a time dependent background it is not meaningful to 
calculate S-matrix elements. The absence of a global time like Killing vector does not allow 
us to relate aymptotic in and out and states and at any rate we mean to calculate the matrix 
element between time evolving states at a given time. The Schwinger-Keldysh |26l . 



27|,|28|,|29| 



formalism addresses these issues by calculating the evolution of matrix elements from given 
initial conditions to some later time. The application of the formalism to de Sitter space is 
discussed in 

We can summarize the results by considering the evolution of the density matrix. If we 
define the evolution operator so that 

i^Uj{r],r]o) =HjUi{r],r]o), (19) 

and we take the full Hamiltonian to be given by 

H = Ho + Q{v- m)Hi. (20) 

We define initial conditions at ?7o where p{rio) = po- Then we can calculate the time- 
dependent expectation value of a general operator, 

^ Tr [Uj{-oo, 0)Uj{0, v)OUiiv, -oo)po] 
^ ^^"^^ [Tr][Uj{-oo,0)Uj{0,-oo)po] ^ ^ 

The numerator of this expression represents a closed time contour (fig. running from the 
infinite past with given initial conditions, enclosing the operator inserted at 1], the infinite 
future and finally running back to the infinite past. 

The calculation is implemented formally by considering fields that live on the forward 
contour and a separate set of fields that live on the backwards contour, labelled respectively 
by '+' and '— '. Time ordering is then implemented so that '— ' fields always appear after 
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countour ^ 



- countour 



I I 

no n 



FIG. 1: The closed time contour that appears in the evolution of operators over a finite time 
interval. Separate copies of the field appear on the forward running and backward running portions 
of the contour. Here the contour is labelled in conformal coordinates. 

'+' fields and that time runs backward on the backward contour. The respective Green's 
functions are given by, 

Q{v - v')Gk{v, v') + 0(V - v)Gk{v', v) (22) 
0(V - v)Gk{v, v) + Q{v - v)GkW, v) (23) 

Gkiv,v') (24) 
Gk{v',v) (25) 

where we notice that G^^ is the usual Feynman propagator. Further, due to the doubling 
of fields, we have '+' and '— ' vertices corresponding to each of the vertices in the original 
interaction lagrangian. Further details of this formal procedure will become apparent when 
we present the calculation of the decay rate. 



III. INTERACTION LAGRANGIAN 

We now describe the model we will use to study inflaton decay |3l which consists of the 
inflaton field coupled to a massless scalar field, a via a three point interaction. We anticipate 
the requirements of the Schwinger-Keldysh formalism and include both '+' and '— ' fields in 
the lagrangian, 

£ = i(9^$+5^$+ + m|<l>+' + d^a+d^a+ + i.Ra+^) 

+ h{t)^ (26) 
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where the negative sign appearing before £($ , o" ) is a result of the time reversal on the 
negatively directed time contour, 

/O p—oo 

r ° (27) 

= - / [Hi{^+,a+)-Hi{^-,a-)]dri. 

J — oo 

The interaction term is apparent, but the magnetic term involving h{t) requires some expla- 
nation. Since we are primarily interested in the decay of $ into cr we need not be concerned 
with the details of the potential for $ as long as the field is rolling slowly with respect to 
the decay rate into a. So instead of including an explicit, but arbitrary, potential we simply 
include a magnetic (or source) term to set initial conditions such that $ is 'slowly rolling'. 
In addition h{t)^ term is required for the renormalization of tadpole diagrams 19|. 

As usual the inflaton field is separated into a homogenous term and scalar fluctuations. 



= 0(t) + where (<l>^) = 0(t). (28) 

This implies the tadpole condition 

(<^±(f,t))=0 (29) 

from which we obtain the equation of motion. In addtion we also enforce (o"(x, t)) = which 
ensures that the a field does not acquire an expectation value. 
Rewriting the action in terms of the shifted fields we find 



S = J d^xV^[Coiv+) + CoicT^-) 



(30) 

+ |(0 + y.+)(a+)2 + /iy.+ 



- [+ -J 

where Cq denotes the lagrangian of the corresponding free theories. Now evaluating the 
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tadpole condition 

= J Sydt^f^i ii) + 3H<j) + ml(f)- ^{a^){t) - h{t) 



.9' 
i- 



dt'yf^K{t,t')(t){t') 



X 



(31) 



47r2 

(((^+(f,r)(^+(y,t)) - ((^+(a-,r)(^-(y,t))) 
+ 0{g') 



where we have defined 



poo 

K{t, t') = / dkk^ {Gk{t', tf - Gk{t, t'f) . (32) 
Jo 

To enforce the tadpole condition we must require that the coefficient of (^{ip~^{x, T)(f^(y, t)) — 
{ip~^{x,T)if~{y,t))) vanish and this yields the equation of motion for 0(i). First we must 
take care of the (a^) term which we simply absorb into h{t) defining 

K{t)=h{t) + ^-{a'){t). (33) 

The inhomogenous term does not enter the expression for the decay rate and so is not of 
concern here. Its main function is to act as a Lagrange multiplier to enforce the initial 
conditions at t = 0, 

m = <i>i (34) 

0(0) = (35) 

(36) 

It turns out that the kernel K{t,t') can be cast in convolution form K(t — t') and so the 
equation can be solved by Laplace transforming (p(t). In general we define 

POO 

/»=£{/(t)}= / dte-^'m. (37) 

The exact form of the kernel depends on the nature of the couphng of the a field to cur- 
vature and its mass. However, generally the kernel integral possesses a logarithmic UV 
divergence which can be absorbed into a mass redefinition. In the Bunch-Davies vacuum 
the renormalized mass is given by 

mlR^ml + ,^{\ne+l) (38) 
10 



The details of Laplace transforming the equation of motion also depend on properties of 



the o field but as an if 
Bunch-Davies vacuum 



ustrative example for the conformally coupled massless case in the 
19^ it is found that 



- + m{s4>{s) - 0,) + m| + -|— E(s)0(s) = K{s) (39) 

S(s) = -si:{ln(l - e-^*) + e-^*}. (40) 

Solving this algebraic equation for 0(s) and taking the inverse Laplace transform the 
decay rate can be identified from resulting solution for 0(t). In the particular example we 
are discussing the decay rate is found to bejl^ 

^^jqmM^^,,,J»ji}i^\ (41) 

167r2 M$ 327rM$ \ 2 ) ^ ' 

where Ph = ^ and to this order 

Ml = ml - ^ (42) 
3H 

s+ = -—+iM^. (43) 

The same general procedure is employed in the following sections. However, the mode 
functions and hence the Green's functions vary for different theories (conformally cou- 
pled/minimally coupled) and a- vacua. This implies that the kernel acquires new features 
and, in the case of a minimally coupled field, infrared divergences. However, after some 
work the decay rate can be extracted in the same fashion from the solution to the equation 
of motion. These results are pressented in the next two sections. 

A. Conformally coupled massless a 

Specializing to the conformally coupled massless case, the calculation for general a pro- 
ceeds in the same manner as detailed above for the Bunch-Davies vacuum, except the ap- 
propriate a dependent Green's functions are employed in computing the diagram. Defining 

= ^ ^ (44) 

a I a* 

= ^ \ ^ (45) 
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we can write the a dependent mass renormalization and self-energy as 

' 167r2 ' 

E(s) = E^{T,e{s) + Do^x{s)) (47) 



mlj^ = ml + E^^r^{\Yie + l) (46) 



where 

is a non-multiphcative correction to the self-energy. It can be seen that the mass renormal- 
ization differs from the euclidean case only multiplicatively and it should be noted that as 
a — oo (the Bunch-Davies limit) the multiplicative factor ^ 1. Although the modifi- 
cation to the self-energy is not purely multiplicative in general, the imaginary part is only 
multiplicatively modified. 



Rewriting the multiplicative factor in a suggestive form 

E^ = l + 2(e-("+"*) - 1)-' (50) 

we can interpret the decay rate as being Bose enhanced with respect to the Bunch-Davies 
rate. It should be noted that the Bunch-Davies rate is unique in that it is the minimum 
rate since > 1. 



B. Minimally coupled massless a 



The mode solution for the minimally coupled scalar in the Bunch-Davies vacuum is given 



by 



-ikr] 



U^{r,)^{i-kv)- 



(51) 



and from this we can construct the general mode functions U^{r]). Then the Green's func- 
tions follow and once again the procedure as outlined above can be invoked. However, new 
complications arise, the first of which is an infrared divergence presumably due to the mass- 
less minimal coupling. The divergence can be regulated by imposing a lower cut-off on the 
loop integral. Further, the resulting equation of motion after Laplace transforming is no 
longer algebraic but rather a first order differential equation of the form. 



Po{s)+X{pi{s) + q{s)-^ 



is) = m 



(52) 
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Solving this equation, order by order in A, 



(Pis) 

Ms) 
Ms) 



Ms) + xMs) + oix^) 

fis) 



Po[s) 

Ms) 

Po{s) 



Pi{s) + q{s 



dint 



ds 



(53) 
(54) 

(55) 



leads to the result 



fis) 



fis) 



(56) 



Po{s) + X[pi{s) + q{s)£\n-^^^^^^ 

This is of the same form as the conformally coupled case, at least to 0{g^). The decay rate, 
once again, acquires a purely multiplicatively enhancement, and we find. 



9' 



327rM$ 



tanh 



nm 



$ J 



(57) 



C. Alpha Time Ordering 

The abscence of divergences in the preceding calculation of the decay rate in an alpha 
vacuum is unusual considering that previous investigations seem to indicate 

that non-renormalizable divergences inevitably appear in alpha vacuum calculations. The 
significance of these divergences is still unclear but in certain cases they can be traced to 
pinched singularities in the propagator. 



As discussed in section piB|l it is possible to define an a-dependent time ordering pre- 
scription that resolves the pinched singularities appearing in the Feynman propagator. A 
related procedure is discussed in j^l where non-local interactions are introduced to absorb 
the divergences that arise while preserving the locality of the theory. In fact, the inclusion 
of a non-local (point/anti-pode) interactions that preserve the locality of the theory in the 



25|. So for this 



sense that commutators vanish outside the light cone was discussed in |23l . 
section we will generalize the interaction lagrangian to include these antipodal interactions. 

Before we present the new interaction lagrangian, it is important to note that the main 
result of the time ordering prescription together with our choice for the arbitrary constants 
Aa = 1 and = is that internal propagators are converted to the Bunch-Davies prop- 
agator and all the a-dependence is restricted to the external legs. This is distinct from a 
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field redefinition in that we hold the interaction term fixed when we change the ordering pre- 
scription. A field redefinition would demand the interaction term undergo the corresponding 

r~i 

transformation. This point is discussed in 25] where it is clarified that a redefinition of the 
kinetic term to make it once again appear local would lead to a non-local interaction term. 
The appearance of euclidean propagators on internal lines indicates that the only diver- 
gences that will appear are the ones encountered in the Bunch-Davies vacuum which can 
be renormalized via the usual mass and (if relevant) wave function and coupling constant 
renornalization. 

We write the new interaction lagrangian that includes antipodal interactions as 

where we take x-p = [rj, x) and Xj^ = {—rj, x) and Ajj is a dimensionless matrix describing 
the relative strengths of the interaction terms. To fix the spurious freedom in the definition 
of A we impose 

A-p-p = 1 (^^) 
A = A^ (59) 

This particular choice of conditions ensures that the earlier result is recovered when all the 
free components of A approach zero. 

Now there are a sum of diagrams to be evaluated resulting in the kernel given by 

k'^dkXijXi'f 

For the massless conformally coupled case, = \, only UV divergences appear in the integral 
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which can be regulated as was done previously with e*'^" — ^ ^ik{a+te) _ ^^lis case we find 



e2 + 4(r] - r]'Y 
Ai{r] + 7]') 



e2 + 4(r/'-r/)2 

+4Ap4App[i(-L_ - _L_) - - - ^^)] 



(60) 



-e " "e2 + 4(r/-V)2 
1 . 1 1 

2^ e - 2i{r] + r]') ^ e - 2i{r] - t]')' ^ 

First of all we should note that the non-local interactions mixing point/antipode fields 
- the X-pji^ terms - give rise to a linear divergence in the kernel. More problematically, all 
terms involving the antipodal field give rise to terms of the form 

fO 



d'q'(t){ri')f{r]',ri) 

in the equation of motion which manifestly violate causality, atleast in our framework with 
adiabatic switching in the infinite past. Boundary conditions of this form have been con- 



sidered in the context of elliptic de Sitter space, see |3l| and references therein. Although 
the theory is local in the sense that communators vanish outside the light cone, we are led 
to conclude that the resulting theory is not causal in our framework since we cannot define 
a useful initial value boundary condition. So it appears that the only meaningful choice for 
the interaction matrix is 

^AA = ^VA = 0. 

With this choice for the interaction the prescription for modifying the propagators recov- 
ers the euclidean result for the decay rate. Further, in 2^ it was noted that euclidean self 
energy is also recovered. This, of course, is not surprising considering that in essence the 
prescription sets all to internal propagators to the euclidean propagator. This suggests that 
proper attention to causality - which motivated the propagator redefinition - removes any 
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vacuum ambiguity in an interacting theory. It provides an additional motivation for calcu- 
lating in the euclidean vacuum, namely, it is the vacuum that yields the natural feynman 
propagator since its definition depends of enforcing causality. However, we can still choose 
to calculate in any vacuum but physical results will be independent of our choice. 

IV. CONCLUSIONS 

In this paper we have considered the decay of the infiaton into light scalars placed in a 
general de Sitter invariant a-vacuum. The result for the commonly considered Bunch-Davies 
vacuum, yields a decay rate that corresponds to decay in Minkowski space at a finite 
temperature, (Sh = 2Ti/H. For a general a- vacuum we find a further enhancement that is 
analogous to Bose enhancement seen in stimulated emission with the enhancement factor 
given in eq.(|5(J|). We find that this is the case for both conformally and minimally coupled 
scalars. This interpretation is compelling considering that when a-vacua are projected on 
the Hilbert space constructed from the Bunch-Davies vacuum they correpond to states with 
constant occupation of all modes which would lead to stimulated decay. Of course, this still 
singles out the Bunch-Davies vacuum as it is relative to the Bunch-Davies rate that the rate 
for general a are enhanced. 

We have also considered a recently proposed prescription for treating the a-vacua that 
seeks to address the proper time-ordering of the Feynman propagator in a general a-vacuum. 
The prescription leads to the internal lines of diagrams being converted to the Bunch-Davies 

, and decay rate. 



propagator independent of a. This leads to a self-energy, caclculated in 
as discussed in this paper, that are independent of a. Further, the only divergences that can 
arise from loops are those that are encountered in the Bunch-Davies vacuum and hence the 
theory is renormalizable for all a. This is an appealing result as it leads to an a independent 
scheme for calculation. However, there is still a dependence in the external legs which will 
appear for example in the two-point function relevant to infiation. This seems to be related 
to the manner of initial conditions and is question that we are still investigating. 
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